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Abstract 

Mechanism of di-pion transitions nS — > n'Snirfo = 3, 2; n' = 2, 1) 
in bottomonium and charmonium is studied with the use of the chiral 
string-breaking Lagrangian allowing for the emission of any number of 
7r(K,rj), and not containing fitting parameters. The transition ampli- 
tude contains two terms, M = a — b, where first term (a) refers to sub- 
sequent one-pion emission: T(nS) — ► ttBB* — > 7rY(ra' 'S)ir and second 
term (b) refers to two-pion emission: Y(nS) — > ttttBB — ► inrY(n'S). 
The one-parameter formula for the di-pion mass distribution is de- 
rived, ^ ~ (phase space) \rj - x\ 2 , where x = J , = M 2 
The parameter ?y dependent on the process is calculated, using SHO 
wave functions and imposing PCAC restrictions (Adler zero) on am- 
plitudes a,b. The resulting di-pion mass distributions are in agreement 
with experimental data. 

1 Introduction 

An enormous amount of experimental data on strong decays of mesons and 
baryons is partly used by theoreticians for comparison in the framework of 
the 3 P model pQ, and its flux-tube modifications [2j. The analysis done in 
[3] confirmed the general validity of the model, whereas in [1] results of other 
forms of decay operators have been also investigated in meson decays, and 
in [5] in baryon decays. On the whole, the phenomenological picture seems 
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to be satisfactory for the 3 Po model with some exclusions discussed in [3] 
and [6]. The extensive study of strong decays of strange quarkonia based on 
the 3 Po model was done in [7J, for a history and review of strong decays see 
Pi- 

Another form of channel coupling Lagrangian was used by the Cornell 
group [JU], and recently exploited in higher charmonia levels [11] . 

Of special interest are the OZI allowed strong decays of heavy quarkonia, 
e.g. 1) ip(n) -> DD, ip(n) -> D*D* etc. and 2) ip(n) -> DD{ixtt), or 3) 
ijj(n) — > DDir(r]). Pionic decays are especially important since they provide 
a fundamental probe of decay mechanism. 

The dipionic decays of heavy quarkonia were first discovered in [T2j , and 
since then attracted a lot of attention, because they display characteristic 
process of a pion pair creation from the nonperturbative gluonic vacuum. 

The resulting dipion spectrum was reconstructed from the general require- 
ments of Lorentz invariance and PCAC [TBI [H] and the multipole expansion 
of gluonic field was used as derived in [TU EES! [16l [T7] . 

In this derivations it was assumed that decaying heavy quarkonium has a 
small radius so that the multipole expansion is applicable. In practice, how- 
ever, the radii of T(2S), T(3S), T(4S'), T(5S), are respectively 0.5,0.7,0.9, 
1.1 fm while these of ip{2S) and -^(3770) are 0.8 fm, [18]. These values are 
larger than the gluonic correlation length of the vacuum, A = 0.2 fm [T9l [20] . 
therefore not the field strength but rather string tension comes into play and 
the internal structure of the decaying state can be important. Therefore one 
should use another and more complete formalism to calculate the transitions. 
As a possible hint in this direction may serve the anomaly in the pionic spec- 
trum. Indeed, while the T(nS) —>■ T(n'S)7nr transitions with n — 2, n' = 1, 
observed by CLEO [2D [221 E3J [21] , show a high mass peak in mr spectrum, 
the T(35) — ► T(1S)tttt transition exhibits a double peak structure. The 
same type of structure was found by BaBar in T(AS) — > T(2S')7r7r [25], while 
both Belle [26] and BaBar see only high mass peak in T(45) — > T(1S)tttt. 
Some theoretical explanations of this anomaly were suggested in literature 
a) the role of final state interaction and a resonance [23 [2EJ [221 ED] , exotic 
T — 7r resonances [23 ED E21 coupled channel effects [Ml E5], relativistic 
corrections [36], S-D mixing [37]. The role of the constant term was studied 
in [38]. 

Another interesting example of the important role of (QQ)tttt channel 
is provided by the X(3872) resonance, which was seen in this channel [39] 
as well as in the channel (DDk) [10], and by recent finding of the F(4260) 
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resonance in the channel J/ipir + 7T~ [41J. 

Recently a detailed analysis of dipion transitions among T(3S f ), T(2S), 
T(1S) states was done by CLEO Collaboration [12], which gives additional 
information on mass and cos 6 distributions, calling for theoretical explana- 
tion. 

From all this set of accumulated data one derives the impression that 
one should construct the theory, where large distances are under control and 
all channels like DD, DD*, DDtt, DDtttt for cc and similar ones for bb states 
should enter on the same ground. Moreover one should try to derive it from 
QCD with minimal number of assumptions and parameters. As a result 
the theory must be applicable to large distances and radii of bound states, 
R > A = 0.2 fm. 

It is the purpose of the present paper to make some progress in this 
direction using the Field Correlator Method (FCM) [13] and background 
perturbation theory [H] to treat nonperturbative (NP) QCD contributions 
together with perturbative ones. 

Physically, the main mechanism is the breaking of string, connecting 
heavy Q and Q, by a light qq pair with a possible emission of Nambu- 
Goldstone mesons. As we shall argue below, there exists a general mechanism 
for the processes of this kind, which is derivable from QCD using FCM, along 
the lines first treated in |45j . 

To this end we shall find the Green's function of heavy Q and Q quarks, 
propagating in the nonperturbative gluonic vacuum, where creation of the 
light qq pair is described by the quark-pion Lagrangian (action) Sqm obtained 
via the chiral bosonization procedure [16] . It is essential, that our derivation 
does not contain any arbitrary parameters, and the only mass parameter 
entering Sqm, M br can be in principle computed via field correlators. The 
resulting structure of the decay amplitude for the process (1) resembles to 
some extent the 3 Po model, and this fact can be used as a step in the 
systematic construction of the theory of strong decays, where 3 Po (or its 
modifications) is participating as an example. 

The paper is organized as follows. In section 2 the general formalism of 
field correlators for the QQ Green's function and for the effective Lagrangian 
of light quarks is introduced. In section 3 the bosonization procedure for 
the effective Lagrangian is described and the effective quark-pion operator is 
written down. In section 4 general relativistic construction of decay ampli- 
tude is given. Section 5 is devoted to the quarkonia decays without pions. 
In section 6 the dipion decays are considered and the expressions of decay 
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amplitude, dipion spectra and total width are given. Section 7 is devoted to 
results of analytic and numerical calculations and comparison to experiment. 
Section 8 contains discussion of results in comparison to other approaches 
and experiment. Section 9 is devoted to summary of results and outlook. 
Acknowledgements are placed in section 10. Four appendices cover neces- 
sary technicalities for relativistic decay amplitude, kinetics and details of 
SHO matrix elements. 



2 Effective quark-pion Lagrangian 

We start here with the standard Euclidean partition function for 3 light 

(mi \ 
m 2 , and ip operators 

m 3 / 

ip 9 , g = 1,2,3 and one heavy flavor with mass tuq and wave operator ^>q 

Z = J DADrD^ 9 D^ Q D^ Q exp[-(S + S 1 + S mt + S Q )). (1) 
We have omitted in flTJ gauge-fixing and ghost terms, and defined 

S = l -J rf 4 x(F^) 2 , Si = -i J d 4 x^ f + m f W f , (2) 

S int = - J d A x$ f gAiP f , S Q = -i j d 4 x^ Q (D + M Q )^! Q . 

To derive the quark-meson Lagrangian one can follow the procedure given 
in [16] . The first step is the integration over gluonic fields DA using the 
generalized contour gauge |47J to express through and the Gaussian 
approximation to cluster expansion as in 

In this way one has 



A fl {x)= I dT llvp {x 1 z)F vp {z) (3) 

J L(x,xq) 

where the contour L(x,xq) connects some arbitrary initial point xq and final 
point x and the measure dT^p is given in 06] and will not be used here. 

The partition function Z does not depend on the choice of contour L, 
and one can introduce the unity operation of integration over some class 
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of contours {L} with the weight Dk{L) In the Gaussian approximation in 
the confining regime the area Sqq inside closed world lines of quarks and 
aintiquarks will depend on the choice of contours {C}, and we shall assume 
that the contour integration Dk(L) results in the fixing of the minimal surface 
Sqq. Actually this minimal surface appears in the area- law asymptotics 
of Wilson loop when all field correlators (and not only Gaussian ones) are 
taken into account, and hence integration Dk(L) restores the action of all 
correlators. 

In what follows we shall be interested in the processes of light quark 
creation, where heavy quarkonium is involved as a background. To this end 
it is convenient to integrate <&/<& iri D^qD^/q with ® in j = (^ Q r in j^ Q ) 
before integrating over DA, which yields 

Z = J DADiPDtPe-( So+ ^ +Si ^G Q Q(x in ,x f ;A) (4) 

where we have defined 

G QQ (x in ,x f ,A) = G^ + G Q i ^ (5) 

with 

G QQ U) = J^ tr { r f S Q( x in, xf, A)T in S Q {xf, x in ; A)}. (6) 

Here GqA refers to the OZI violating decay and production processes, which 
will not be discussed in the present paper, and this term will be omitted in 
what follows. 

The connected part of Gqq can be written using the Fock-Feynmann- 
Schwinger (FFS) path-integral representation [IS] 

G { QQ n \x f ,x in ,A) = J dp(C Q )W(C Q ;A) (7) 

where the Wilson loop operator is defined as 

W(C Q ,A) = ^-tr c P A (exp(ig I A^dz,)) (8) 

iV c JCq 

and the dp{Co) includes integration over paths C\,C2 of heavy quark and 
antiquark dp{Co) ~ DzDz, so that Cq — G\ + CV 



5 



The next step is the integration over all gluonic fields in (j3J) which re- 
duces to the expression f DAe~ Sint W(CQ, A). Here enters linearly in the 
exponent and one can use the cluster expansion theorem and keep only the 
quadratic terms in A^ (and due to ([3]) also quadratic in F^) - this approxi- 
mation which is usually called the Gaussian one - is justified by the Casimir 
scaling accurate to few percent [15] . 

The result of this integration can be written as (see [5Uj for details of 
derivation) 

DAe {SmtiA)+SoiA)) W{C Q , A) = (e~ Smt{A) W{C Q ,A)) 



= J DK(L)e~ s *ff^W(C Q ,L). (9) 

Here S e ff(L) is the same as in the case of light quarks without heavy quarko- 
nia and W(Cq, L) are expressed through the same field correlator. 



S e ff(L) = -- I dxd y$bl^a(x)$a'(y)lwfo(y)(8aa'8w ~ J^^ab^a'b')J , 



(10) 

w{c q ,l) = exp ( - ^y dx » J d yuJ^u(x,y)j (11) 



'Cq JCq 

where J^ v is 

J^(x,y)= dT, Xp ( Xj z) / dT^ v (y,z')D Xp ^ v (z,z'), (12) 

J L(x,xq) J L(y,xo) 

with the gauge-invariant field correlator 

g 2 

Dx P ,^(z,z') = — (trF Xp (z,x )F iv (z,x )), (13) 
expressed through the parallel transported field strength, e.g. 

F piy (z,x ) = $(x ,z)F IJi1/ (z)$(z,Xq),$(z,x) = Pexpig / A p du p . (14) 

J X 

For illustration purpose, when only heavy quarks in Gqq are present, i.e. 
S e ff = in Qj, the point xq can be chosen inside Cq and contours L(x,x') 
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taken as straight lines, L — > Lps like in the Fock-Schwinger gauge) and W 
acquires the standard form: 

W{C Q ,L FS ) = exp(-~ I da Xp {u) I da iri (v)D Xp ^ v (u,v)) (15) 

where da\ p (u) is a surface element at the point u, S m i n is the minimal area 
surface inside Cq, when confining part D is present in D\ P £ V , 

D\ p £ V (u, v) = (Sx(S pv — 5\ v 5 P £)D{u — v ) + total derivative (16) 

then ( 1T51) reduces for large contours to the area law form, W(Cq, Lps) — 
exp(-aS min ). 

Now consider the case when both terms in are present. To the lowest 
order in S e ff(L) and integrating over light quarks, one obtains in Qj a factor 
(SeffiL))^ = -if d A xd 4 y{-i p S{x,y)'-{ v S{y,x))J pv (S(x,y) is the light quark 
propagator) similar to Gqq, Eq.([6]), which produces as in ©,([8]) a new Wil- 
son loop of light quarks. The situation here is similar to the calculation of 
the vacuum average of the product of two Wilson loops, which was done in 
[5T] . Indeed both exp(— S int (A)) and W are linear in A p (hence also in F pv ) 
and the general formalism used in [51] is applicable also to calculate (J9]). In 
case of the connected average of two Wilson loops, 

x(Ci,C 2 ) = (W(C h A)W(C 2 ,A)) - {W{C X ,A))(W{C 2 ,A)) (17) 

the resulting x{C±, C 2 ) for the case of two large contours Cj, C 2 in one plane 
with opposite orientation was obtained in [5T] for 5™ m < S*™ 1 " 

X (C 1 , C 2 ) = exp(-a(5 1 min + Sf n )) exp(2aS 2 min )+ 
+ (l - J*) exp (-^) - l} = i [exp(^(r - 5 2 min )) + O 



7V|-1 
(1 , 

The net outcome of the analysis in [SI] for a general configuration of two 
oppositely oriented contours C±, C 2 is that for small distance h between the 
surfaces Sf 1111 and S 1 ™ 111 , h < h crit one can write 

x(C 1 ,C 2 )^-^exp(- C r5 min (l,2)) (19) 

c 

where 5 imn (l, 2) is the minimal surface with the boundaries at C\ and C 2 , 
and the critical distance h cr u is defined by the condition which can be ap- 
proximately written as such distance h for which S min (l,2) = Sf 1111 + S 1 ™ 111 . 
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For h > h cr it the value of x( c i> £2) is defined by perturbative exchange of two 
gluons propagating in the surface S?>^ m (l, 2) (see [5T] for more details and 
discussion) - the two-gluon glueball exchange between contours C\,Ci. 

Coming back to the calculation of (Q one can notice that the same answer 
f fT8|) for two Wilson loops can be obtained by the proper choice of the set of 
contours {L} which minimize the common surface S(l,2), namely 

x(Ci,C 2 ) = -L(exp(- < 7S , L (l > 2))) ininL . (20) 

Indeed one can show that the contribution of the kernel J^ v to the product 
of two oppositely oriented Wilson loops on one plane Wi(Cq, L)W / 2 (Cq, L) 
vanishes inside the smaller loop W%, which yields the same answer as in ffT9]) . 

In a similar way the average in ([9]) can be written as 

(e- s in*wW(C Q ,A)) = fD^L)e~ s ^ff^W(C Q ,L) = 

= ^e- s ^ L ^W(C Q ,L ) W 

where Lq is the set of contours which ensures the minimal area between the 
quark trajectories generated by S e ff (those will be exemplified below) and 
trajectories of heavy quarks. 

As a result of a light quark pair production in presence of the heavy quark 
loop W(Cq, A) one thus obtains after averaging over vacuum fields the same 
loop covered with the film, but with a hole due to light quark loop. This 
situation is depicted in Fig 1. Now we turn to pion creation in the same 



system. 
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3 Light quarks and pions in the heavy quarko- 
nium 



Having in mind the main result of the previous section, Eq. (|2ip . one can 
proceed as in the case of light quarks without heavy quarkonium, i.e. as was 
done in [46J, correcting in the final results for the holes, made in the world 
sheet of the heavy quarkonium by light quark loops. In this way we first 
obtain as in [46] due to the bosonization of the four-quark action S e ff(L), 
Eq. ([TO]). 



Z = j Di;Di,e~^ +s ^DM s D(t )a d P {C Q )W{C Q ,L) (22) 
where the quark-meson effective action Sqm is 

S QM = - J d 4 xd 4 y$f(x)iM s (x,y)Uf9(x,y)r(y)- m) 
-2iV / (J Mi/ (z,y))- 1 M s 2 ( 2 :,y)], U = exp(t l5 t a Mx, y)). 

Integrating out the light quarks one obtains the Effective Chiral La- 
grangian (ECL) as in 



with 



where 



Z = J DM s D<p a dp(C Q )W(C Q ,L)e- SEC \ (24) 
S ECL = 2N f J d 4 xd 4 yJ^M*(x, y) - W (</>), (25) 



W(<f>) =N c trln[i(d + rh + M s (x,y)U)} (26) 



The integration over DM s D<p a is done in a standard way using the stationary 
point equations in ff24|) which yields 



0i°) = O, Mi°\x,y) = ^-J^(x,y)Tr(S(x,y)) (27) 

where S(x,y) = S^x, y)|^= , and 

S+{x,y) = (xKtd + im + iM^U)- 1 ^}. (28) 
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In what follows we shall be interested in heavy quarkonia decays via the 
OZI allowed processes with emission of zero, one and two (in principle, more) 
light mesons. 

The latter are present in the factor U in fl26|) . which can be rewritten as 



U = exp(i 75 a t a ) 



exp O'^^J^j ' V?aA a = v 7 ^ 



/ -4= + — 7T+ K + 

(29) 



and = 93 MeV. 

Moreover the asymptotic solution of stationary point equations (I2T1) . per- 
formed in [52] yields (|x — xo| A) 

Mf»(x,!/) «<r|x-xo(L)|5< 4 >(x-y) = M(-x)5 {i \x-y) = M(x)5 (4) (x - y) 

(30) 

where x (L) is defined by the choice of set of contours L, and as we argued 
above, the minimal set of contours corresponds to the minimal area of the 
surface between heavy quark trajectories with the hole made by light quark 
loops made of S^x.y). One can easily understand that this situation is 
realized when x (L) for a given light quark is chosen at the position of the 
heavy antiquark, and correspondingly for the light antiquark. 

One important conclusion from (|30|) is that for large distances R ^> A, 
one can use only local quantities for and (p a ,(p a (x,y) — > (p a (x). 

As the next step we expand W{4>) in powers of the Nambu-Goto fields 4> a 
and consider separately terms of the zeroth, first and second power in <fr a in 
the partion function Z. 

One has 



W{(f>) = NjrlnlS- 1 - M T5 - |M0 2 ] = 

w (<p) + Wi(<!>) + w 2 {(f>) + ...,4> = t a <p a = ^ 

Jit 



(31) 



W 2 {(j>) = ~ ytr(^M0 2 + S.\lo~:,S~:,.\lo). (32) 

In an analogous way one obtains terms of higher order in <f), which corre- 
spond to the decay processes of heavy quarkonia with emission of three and 
more NG mesons. 

At this point one can realize, that the interaction between heavy quark 
loop and light quark (and pions) is mediated by M s (x,y), which is string- 
like for large distances (large loops), as it is in physical situation. For small 
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loops, one has instead in (TT21 (M s is proportional to J^ v ) gluonic condensate 
(F^ u (x)), as it is in [TUGD]. Below we shall discuss the large loop situation, 
and to this end we need to detalize the hadronic states in the big heavy quark 
loop with the hole made by light quark loop. 

Consider now the two-pion term W2(4>) in the background of the heavy 
quark loop W(C Q ,L), as in f[2"l"l) . 

The amplitude is proportional to 

(W 2 (<j>)) Q = I dp(C q )W 2 (<P)W(C Q ,L) (33) 



One can rewrite W 2 (4>) using [53] as follows 

W2i<P) = \\ $^^ Mh)N{ku k2)Mk2) (34) 

with 

N{hM) = y {/ dxe l ^ +k ^ x tr(AM s ) xx + 
+ J d 4 xd 4 ye iklX+ik2y tr(A(x,y)M s (y)A(y,x)M s (x))^ (35) 

and 

A = (d + m + M s )-\ A= (d-m- M s )-\ (36) 

Both terms in (1351) are depicted in Fig.2, the first-tadpole term in Fig. 
2a, and the second-two-point quark loop in Fig. 2b. 




Fig 2 (a) Fi g 2 ( b ) 



The resulting dipion spectrum is defined by N(k, k') averaged with the 
heavy quark Green's function (only one part of it - the Wilson loop W(Cq, L) 
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is present in f)33p for simplicity of discussion) , so that one writes for the dipion 
spectrum with the phase space factor 

dV-KTv ~ \ (N(k 1 ,k 2 ))Q\ 2 d^(X7T7r)- (37) 

At this point it is convenient to discuss the PCAC limit of our expressions 
( |33|) . fl35|) . e.g. the so-called Adler zero problem. In terms of N(k,k') this 
amounts to the requirement of vanishing of (N(0,0))o in the chiral limit, 
m q = m — > 0. 

In the vacuum average case, when no heavy quarks are present, this 
property was proven in [53] . 
Indeed, one can write 

(iV(0,0))o = i(tr{(AM s )o + J d*ze lkz A(0, z)M s (z)A(z, 0)M S (0)}) = 

= -{tr(AM s A(d-m))) = -mtrA = ~ W {tnM>) = (38) 

We have used here identical transformation for the quark loop: 

(tr(AM s )) = (tr(M s AA^A)) = (tr(M sl5 A l5 (M s + m)A)} 

exploiting the symmetry x M — > — of the vacuum. 

Thus in pure vacuum one obtains expansion of N(k, k') = 0(m)+0(k 2 , k 2 , kk') 
Note the relative negative sign of two contributions in Fig. 2a and Fig. 2b 
However, in the realistic situation of the reaction X(n) — > X(n')inT the low- 
est order terms 0{ki fJ ,k2 V ) from the expansion of N(k±, k 2 ) are multiplied by 
(among others) momenta of heavy quarks, so that a typical term in the pro- 
cess amplitude would be , where V is the momentum of X(n). In 
the cm. system it amounts to amplitude ~ which according to Eq. 
(1A2.6P is ~ ( Ah p ; where AM is the mass difference of X(n) and X{n'). 
This factor does not provide any noticeable damping of tttt spectrum near 
threshold. Therefore we do not have any reason to believe that the Adler 
zero argument alone can play any role in the building up the nir spectrum. 
An additional support of this statement comes from experiment, where spec- 
tra in T(3S) -> T(15)vr7T, T(2S) -> T(lS)vr7r and T(3S) -> T(2S)tttv do 
not have a universal damping near irn threshold [12], but are quite different. 
At the same time the spectrum of ip{2S) — > J /ipilS^ir [51] behaves like 
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~ (g 2 — 4m 2 ) 2 , as suggested in [53]. This behaviour is not deducible from 
Adler zero arguments. 

However, the phenomenon of two terms with different signs in Fig. 2a 
and 2b, fully compensating each other in the chiral limit, is vitally impor- 
tant for our final results, where these two amplitudes will have destructive 
interference. Moreover, in section 7 we show that applying the Adler zero re- 
quirement to our final equation, one predict the amplitude of im production 
in good agreement with experiment. 



4 Derivation of the decay transition 

The original idea of the mechanism discussed in this section was given in 
[4"5| . and was called there the Chiral Decay Mechanism (CDM). 

To generate the string-breaking vertex one can use the quark-meson La- 
gragian (123]) . which for the mesonless decays can be written as the string- 
breaking vertex operator 

Sqm -> S S Br = -i J d A x J d 4 y^(x)M s (x,y)U(x,y)ip(y) -> . , 
-»• -ij d 4 x^(x)M br U(x)ij{x). 1 ' 

We have introduced in f]39l) the local limit of the mass operator M s (x,y), 
which we call the vertex mass operator M br . It differs from the usual mass 
operator M a (x) = a\r\ placed at the end of the string and incorporating 
confinement, as exemplified in (1301) . because M br is positioned somewhere in 
the middle of the string and at the beginning (or the end point) of the light 
quark loop. Therefore the calculation of M br needs a special care, and was 
originally done approximately in [35] . One should stress that M& r is placed 
at the pion emission vertex of the quark-loop hole in the film of QQ string, 
a similar important role is playing by the vertex mass M(0) in the current 
correlator diagram in [53] and the value M(0) ~ 150 MeV for u, d quarks 
was used to calculate f n and (qq). 

However dynamics in M br and M(0) is different, since the former defines 
the decay process of a long string, while the latter refers to the amplitude 
of qq meeting together at one point, with the (short) string between q and 
q. Therefore M(0) ~ aX, while M br is expected to be of the order of (m + 
D) ~ m q + Co. In what follows we shall consider M br as a free and universal 
parameter, and we shall fix it once from the decay ^(3770) — > DD and use 
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for all dipion decays of T and if) families. As we shall see, this strategy gives 
a good result for the absolute values T^, at least for T(2S) — > T(15')7r7r. 

To proceed we shall define the problem of quarkonia decay with emission 
of a light quark pair and any number of NG mesons. To work in a relativistic 
invariant way we need to define the vertices (currents) for creation of heavy 
quark pair in a given state at point 1, Jqq = ^/qTi^/q, with Ti = j if 75, 757^ 
etc. and the vertex for the creation at the point a; of a light quark pair 
plus possibly some number of NG mesons is given in (|39|) . T x = Mb r U(x). 
It is clear, that the general outcome of this creation will be production of 
two heavy-light mesons at points 2 and 3 with vertices ipQT2ip q and ip q T 3 ^)Q, 
see Fig. 3. The connected Green's function for such production, Gu3x, is 
expressed via quark Green's functions SQ( q )(i, k) and has the form 

G 1 23 X = (0\J 1 J 2 J3J(X)\0) = 

= (tr(r 1 S Q (l,2)T 2 S q (2,x)T x S q (x,3)r 3 S q (3,l))). (40) 

2 




As it is shown in Appendix 1 and 2, one can use for Sq^) the Fock- 
Feynman-Schwinger Representation (FFSR), see [56] for review), where all 
Dirac bispinor structure is factorized in the first approximation (neglecting 
spin-dependent forces as compared to confinement potential). Thus one in- 
troduces the factors Z (see Eq. ( lAl.lip in Appendix 1), and the rest part of 
G\23x does not contain spin factors, but rather overlap integrals of the corre- 
sponding scalar wave functions, when the spectral representation of Green's 
functions is used, see Eq. ( 1A2.21|) in Appendix 2. 



The remaining problem is now to go over from the current correlators, 
i.e. 4 point Green's function Gu3x, to hadron-hadron correlators, which is 
done by amputating the point-to-hadron pieces, proportional to the decay 
constants of given currents, as shown in Appendix 2. 



14 



The corresponding "point-to-hadron" matrix element is given by (see 
(IA2.11I) for details) 

(0\jr\n,P) = er^-^f^, £r = 1 for scalars, ey — for vectors (41) 

where decay constants are computed through solution tp n (r) of relativis- 
tic Hamiltonian, as in [57] . 

i f (n).2_ 6|y n (o)i 2 y r u , 

l/r I - M , 



and 



tr^K-A)!^-^)) 

y rir 2 = — — z , (4dJ 

4cJi,a;2 



i.e. the same as the factor Z for the two-point function, = Zi, see 

Appendices 1,2 for details. 

Now dividing Gu3 X by matrix elements ( 14T|) and taking Fourier transform 



in the coordinates ofjpoints 1,2,3, one obtains the hadron decay amplitude, 
as derived in ( 1A2.23 



S Pl |G|n 2 P 2 , n 3 P 3 ) ee n/ j / | Gl pf 4X) r i v'r 3 n = ( 44 ) 

U{0\ji\niPi) exp(- E n ii) 



(27r) 4 5( 4 )(P 1 -P 2 -P 3 )^ 1 n 2 n 3 (p) 



and defining y 12 3 = z } 23 * , one can write 

J ni n 2 n 3 (P) = 2/123 ^ d 3 (v - u)rf 3 (x - u)e* Pr ?/'* l (u - v)^„ 2 (u - x)^„ 3 (x - v) 

(45) 

where ip n . are solutions of the corresponding Hamiltonians for heavy-heavy 
(ni) and heavy- light quarkonia, see Appendix 1,2, and Eq. ( 1 A 1.141) for de- 
tails. 



1 a. similar procedure is applied in lattice simulations, e.g. in case of semileptonic decay 
formfactors, see [55] 



15 



To proceed we introduce QQ Green's functions, in the energy representa- 
tion with the cm. momentum equal to zero, having spectral representation 
(without light quarks) 

G QQ {^\E)=Y, e_ Q e (46) 

n 

and for the QQ Green's function with insertion of the light quark loop one 
has 

a* 2 . E) _ v ^>>^W<> + (2) , 47) 



where 



(E — E n )(E — E m ) 



^)=7y^ p E £ _ Wp) , (48) 



Here 7 = and factors Z are computed in the Appendix 1 together with 
details of derivation and definition of einbein masses u^, fi^. 

It is easy to derive from fj4*8|) an expression for the (complex ) shift of the 
n-th energy level, which has the same form as in the general multichannel 
theory (see [59] for a review, and [60], [61] for a recent application to the shift 
of D*, B* masses). 

E-E n = -w nn {E) (49) 

and the width in the channel k is obtained as (for two equal mass final 
heavy-light mesons) 

r « = 1 PJ ^T J ^pl^,n 3 (p)| 2 (50) 

where M, pk is the reduced mass and relative momentum in the channel k. 

Expressions (1491) . (I50p are sufficient to obtain shift and width of any state 
due to a light pair creation, provided the dynamics of heavy-light mesons 
is known (wave functions and einbein masses (average energies)). As was 
pointed out before, additional NG meson emission is described by simply 
multiplying integrands of matrix elements J nn2m with terms of expansion 

U(x) = exp (z75%M = 1 + ils^r^ + <-Pa = f^f , adding meson 

\ Iw / Jw \/2u;(k)V3 

energy in the denominator and modifying r in fj45l) . In the next section 
we discuss as a clarifying example the decay without meson emission , e.g. 
^(3770) — > DD and in section 6 emission of two pions. 
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5 Heavy quarkonia decays without pion emis- 
sion 

In this section we calculate explicitly the operator W(E) in (|4"T1) . (JUJ), which 
reduces to the calculation of the transition matrix element J nin2 n3 Q45p ■ As 
a practical example we have in mind the decay T(45) — > BB, or ^(3770) — > 
however our results will be applicable to other two-body decays of this 
sort. We start with the calculation of the factor Zi 23x , given in (1431) and 
Appendix 1, where we should insert Tq = ji,TQg = 75 = T q Q,T x = 1. The 
resulting trace is easily computed and can be written as (Zi 23x = Z) 

Z = tr L [%(m Q + fi Q 7 4 - ip Q )l^{m q - u q >y 4 + ip q ) x 

x (m, + ^74 - W q )ls{mQ - fig7 4 + zp<g)]. (51) 

Here notations are p^ = Yl^iPkl 1 - For the or DD decays in the 
cm. system one can simplify itlq = tuq, Qq = Qq, and one finally has 

Z = 4{i2m Q ttuj(p qi - p#) + iu 2 m Q (p Qi - pgj 

+ im(p Qi (p q Pq) + PqiiPQPq) ~ PqiiPQPq) ~ PQi(PqPq) ~ Pqi(PQPq) + PqiiPQPQ))} 

(52) 

One can define relative momenta in the Qq and Qq systems so that in 
the total cm. system (see Appendix 3 for details) 

PQ ~Pq PQ - Pq ,_„s. 

2 = Qi> 2 = q2 ' Pq P « = P ' p Q + p i = ~ p ' ( 53 ) 

and finally 

, P , P P P fKA\ 

Pq = qi + g ' P? = _c *i + 2 ' p = ^2 - 2 ' P? = _c *2 - 2 • ( 54 ) 
The normalized factor Z = „ — ^7= — r is for large mass mo 
- ^ 4 • 2m Q ittu(q li - q 2i - p») _ im Q (q u - q 2i - pi) 



The transition matrix element </nin 2 n 3 (p) = </(p) is 



(55) 



J(p) = y Ze ipr d 3 (u - v)d 3 (x - u)^+(u - v)^ n2 (u - x)^ n3 (x - v) (56) 
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where r = n( - u t ^ = c(u — v). 

It is convenient to introduce Fourier transform of wave functions, 

*»(») = J^fT^), ^n 2 (u-x) = /^ qi(U " X) ^^ 2 ( qi ) (57) 
As a result one obtains for J(p) 

yi23^p^n(cp + qi)^n 2 (qi)^n 3 (qi) (58) 

and Z simplifies, 

Z(q =q . 2 ,p) 2{ ^_^H +0 (V\ m 



As one can see from (I5U1) , Z is proportional to which signals the P wave 
of relative DD or SB motion, as it should be. 

We turn now to the calculation of the factor 2/123 in (Hoi) , which can be 
written (since M br is factored out in T x and Z, Eq. (I5TI) ). 

7/123 = , (60) 

ZIZ2Z3 



with Z given in (1591) and computed in [57J, 

tr(rj(mi - zpi)ri(m 2 + ip 2 )) 



mim 2 + UJ1UJ2 + Vi(P ) 



(61) 



o;iu;2 

Here r\i — +| for the vector and -L for the pseudoscalar case. Taking 
into account, that with the reasonable accuracy Cdi = (p 2 ) + mf, one can 
compute for the vector (i[)(nS),T(nS)) the scalar (D,B) cases, 



n > *~ a 2 



Z D - p; , F^f ( 62 ) 
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where we have denoted = Q D in D,w e u q , q = u,d and approximately 
(see Tables in Appendix 1) ft D = « 1.5 GeV for m c = 1.4 GeV. 
As result y 123 (ifj -> « Z ■ 1.1 = l.l ( ~^ me . 

To compare the decay probability of ?/>(3770) one can use (J5D]) to write 
(we use (p 2 ) = |(p 2 )) 

One can see the correct p 3 behaviour of the p-wave decay, and the 
f actor for the string breaking effect, as it should be. 

At this point one should use realistic wave functions of if), D to compute 
J(p). Here and below in the paper we exploit for semiquantitative estimates 

the SHO wave functions, = ^i(p) ex P ( — 2^f)> & ven i n Appendix 4 

together with matrix elements, where oscillator parameters are fitted to re- 
produce the known r.m.s. radii of states, (r 2 )^ = (0.8 fm) 2 , (r 2 ) D = (0.57 
fm) 2 162], which yields /fy = ft = 0.467 GeV, ft, = ft = 0.42 GeV and J{p) 
takes the form 

,(-) ( p) = ( (f ) 3/2 („ - e-* (64) 

with A = 2/3 2 + /3| , -y = &&., c=^. 

For the reaction ^(3770) -> one has p = 0.276 GeV and j( s/ra )(p) = 
5.09(GeV)- 3 / 2 . As the result one obtains 

MbrV ( J 



r = 19 MeV — ^ . (65) 

\2luJ V5.09GeV- 3/2 y 1 ' 

Since in PDG r e2 , p = 23.6 MeV one expects in the SHO approximation, 
that 

M br = 2u^l GeV. (66) 

This estimate will be used below for dipionic transitions. (Note, however, 
that our treatment of ?/>(3770) as the 25 cc state is a crude approximation, 
?/'(3770) being mixture of 3 -Di, 3 5*i states, therefore a mixing coefficient 
should appear as a factor of J. In any case our expression (|65p is rather 
qualitative, giving an order of magnitude estimate of the decay process). 
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Exact calculation is sensitive to the details of the wave function and will 
be given in a separate publication. One can say, that calculation of T is a 
good check of the form of wave functions, which can be used in addition to 
radiative transitions and lepton decay width calculations, checking it with 
experiment. As it is, we can now turn to the main subject of paper, having 
a formalism without fitting parameters. 



6 Heavy quarkonia decays with emission of 
two NG mesons 

We shall consider in this section the string breaking in heavy quarkonia with 
emission of two NG mesons. Here we come back to Eq. (l23]) for the string 
breaking action and using it one has 

Ssbv = ~i I d T1 xip f (x)M br U f9 iP 9 (x) (67) 



where we have made explicit flavor indices, and U is given in (12"9"|) . In the 
previous section we have considered decay without NG mesons, and corre- 
spondingly replaced U^ 9 by 5f g and Mf\x) by the vertex mass M br . Here 
we take into account higher terms of expansion of U, namely 

U = exp ^*75-^— j = 1 + *75-^ — + ... (68) 

We shall disregard in the first approximation interaction of emitted NG 
mesons with other particles (taking FSI into account as the next step) and 
write 



¥a{x) = -7==, uj a = Jk 2 + m 2 a . (69) 

V V3 

As before one can define the QQ Green's function containing one light quark 
loop with NG mesons emitted from loop vertices at points x and y, which is 
again has the same form as in ( 1471) . e.g. for two pion emission one has 



„^(l)«)i* m * , (S)VS )+ (2) 
Ggg-a^-S %_ En) l E '_iy - (TO) 
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But Wnm(E) now consists of several terms, corresponding to the diagrams 
depicted in Figs. 4 and 5, namely 




(71) 



Fig 4 Fi § 5 



Note, that in (ITT]) indices n«n^ and n" refer in general to three different 
sets of intermediate states; also S(ki,k 2 ) = \/k^ + ml + ^/k^ + m\ . 

Now J n n 2 n 3 (p) and 7 in (17X1) are the same as before and defined by (F4oT) 
with Z ns defined in ( lAl.lip and r x = 1, the dash sign refers to flavor indices). 
For Jnn 2 n 3 (p, ki) one should take into account ([58]) with extra factors of 
^j^-(fa and pion plane wave (ESI). 

As a result one can write 

r _ «pr+iki(x-R) 
J { X n3 {vM) = J Z^ 2ns d 3 (u-v)ci 3 (x-u)vI/+( U -v)^ 2 (u-x)^ 

(72) 

where Z^j} is now also a trace over flavor indices, 

^q^^q-f-(VaK)-Ki (73) 
Jit 
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and the notation (ip a ^a)n implies the numerical matrix, obtained from (j29p for 

/0 1 0\ 

a given pion, e.g. for tc + it is \f2 I , while for ir° it is ((pa^a)^ = 

\ / 

\ 

—1 , so for SU(2) isospin group one has ((f a ^a)n = = tt + t_ + 
/ 

7t~t + + 7r°r 3 , while Tj Pauli matrices. For isospin conserving decays also 
one of vertices Tq^, TQ q or both, should have nontrivial flavor structure, 
otherwise the flavor trace will be nonzero only due to quark mass matrix 
/ m u \ 

m = I rrid , which enters in the light quark Green's function 
V m s ) 

(|28l) . This is what happens in the decays of the type T(35) — > T(lS)r}(ir), 
and will be considered in a separate publication. 

To complete the flavor trace in Z n k one should specify the isospin structure 

( B+ 

of intermediate heavy-light mesons, e.g. for the isospinors B = 

V -°o 

and B = ( ^_ j one actually has final states in the combination K% 

ni(BTiB) and this factor appears in Z n k for given final state of BB (or DD 

'{< 



is the case of 



(2) 

We now come to the vertex with emission of two pions, Jnn 2 n 3 {p, kj., k 2 ), 
which is due to the presence of the factor U in the vertex, specifically of its 
quadratic term, 

Maf {^r^n) 1 

9f 2 = ^f2 = ^TT^i 77 ! 77 ? + ie ikl^lik^2kn)- (74J 

^Jtt ^Jtt ^/tt 

One can notice, that in the case of one-pion vertices in the first two terms 
in (ffll . the resulting isospin structure is such that pions appear in the total 
isospin 1 = state, unless some charge (isospin) filtering of intermediate 
states is done, indeed 

K^Kp n) = Tvu ^(BTiB)(BT k B)ir 2 k = tr f i {-KuTi^kTk) = 2ttitt 2 . (75) 

B,B 



In contrast to that, the Jnn 2 n 3 term contains the 1=1 term, which 



r(2) 
nn: 

of course vanishes for the transition between 1 = states, as in T(ns 
T(n's)TTiT, but can be nonzero for BB final state. 
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As the result, one can write Jnn 2 n 3 as 

- (7nr) e ^(k 1+ k 2 )x rf3(u _ v)d 3 (x _ u)l ^+ (u _ v) ^ 2 ( u _ x) ^ (x _ v) 

V2^(k 1 )y 3 2o; 7r (k 2 )V r 3 

(76) 

and Znn^nz is obtained from the general formula (1A1.11I) substituting — ► 



_ 7Ti7r2 



. In the same way one can consider emission of more pions. 



2/1 

Finally one can define the decay probability for the process (QQ) n — > 
{QQ)n>TT"n, which is obtained from the amplitude w^,\e) by the standard 
rules 

J ft \ l l\ \ I M/ tti\ k ^3<^ 3 k2 

- (n )mr) = K n ,>(E)\ 2 ^ ^ x 

x27r<J(#(k 1 ,k 2 ) + ( 77 ) 

It is easy to see, that the factors V3 cancel in ( 1771) and diu on the l.h.s. of 
(177)1 has the dimension of mass; integrating over IT<i 3 kj on the r.h.s. of (J77I) 
one obtains the total width in this channel 

Tt n ' ] = f dw((n) -> (n>7r). (78) 

Thus one obtains in (jTTjl . ([75]) the absolute rate of the process and fixing the 
dipion mass 

q 2 = M 2 W = (h + k 2 f = 2ml + ^1^2 ~ 2k x k 2 (79) 

one defines the dipion spectrum measured in experiment. 

To proceed with the analysis we consider Jnn 2 n 3 (p, ki) given by (172"j) and 

— (it) 

to calculate the Zkn 2 n 3 , one should specify the intermediate state n 2 n3. It is 
easy to see, that the closest threshold is given by the S'-state of (BB* + BB*) 
plus a pion. Denoting by index k the spin direction of B* and by index i 
the spin direction of incident (QQ) state, one easily calculates zj^ in the 
approximation, when all momenta are small compared to the heavy mass 
rriQ,0(^). The result is (for the intermediate state BB*) 

Z lk -x— J^S*—xp- ■ ( 8 °) 
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Now introducing Fourier transformed wave functions as in (j58p . one arrives 
at the final expression 



f d 3 q\ ki 

(81) 

Comparison with (1551) shows that Z in (THTT) corresponds to the S-wave, 
and dependence on kx appears in the arguments of wave functions, which 
influences the resulting spectrum of nn, however the resulting dependence 
appears to be rather moderate. 

We now turn to the term Jnn 2 n 3 (p, k 1; k 2 ) in (17T1) . the general form of it 
is given in (ITS!) . 

First of all, we calculate ZnnXzi and define the intermediate state BB in 
the P wave - the same as in the decay without pions. Hence the final form 
of Z is 



^4if)(-sK--> 



-TS (A -- p ' ) - Kski+k2 ' (82) 

Here we have used Appendix 3 to express momenta in terms of K, p, and 
Finally J (2) can be written as 



(27r)V2^(fci)^ 3 2^(fc 2 )y 3 



rf 3 gi* n (cp + q x + ^)^„ 2 (qi)^n 3 (qi - K) 

(83) 

Note, that in (1711) this matrix element is multiplied with another one, where 
pions are not emitted, and it has the form (158|) . (159|) with \& n / corresponding 
to the final state of heavy quarkonium. 

Equations (15^1) . (18~T1) and (1831 give all the necessary elements for calcula- 
tion of 7T7T spectrum and total width of two-pion transitions between states 
of heavy quarkonia , provided wave functions of initial and final QQ states 
and heavy-light mesons Qq and Qq are known. However, the resulting mul- 
tidimensional integrals are not easy to evaluate, and we resort at this point 
again to the SHO wave functions, as we did it in previous section. For the 
case when all wave functions are taken as in oscillator potential, the SHO 
integrals for J, , can be calculated in analytic form. 
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Indeed, using Fourier transforms of oscillator wave functions as in (I58p . 
(JEU), dHSD and writing 

*n(q,P) =V n (q)e~^ (84) 
where VM = {^ff\ P M = ft (l - , 

t>M = (¥) 3/2 4i ((f) 4 ~ 5 (if + fj ■ one easily inte g rates over 

d 3 qi in (EED , dHU) , dH3j) with the result 



P 2 

J(p) = Z(p)e ^?/„(p) 



(85) 



D 2 

and J n (p) is a polinomial in with coefficients given in (A4.2), (A4.3), with 



jW(p,k) = lP ' j e ^? ^/ w (p). 



(86) 



Here ZW(p,k) is given in ( ISTfl) . Finally the 27r matrix element is 

P 2 K 2 



J (2) (P,ki,k 2 



v /2^(A; 1 )V r 32^(A; 2 )y3 



4(p) 



(87) 



and Z (7r7r) is given in (152)1 . 

Insertion of (1551) . (1561) . ( 1571) into (17l1) yields the form (i,k stand for po- 
larizations of initial and final QQ states) 



nn'ik\ I 



7 



^2uo^k 1 )V- i 2uo^k 2 )V z f, 



m. 



x 



k'+k 2 

+ ^ 2 \ 2 (^AJ^^AJ^e^^T 



2ft 2 



£ - E'(p) - £(kx) 



+ (1-2) 



K 2 

4^f 



+ 



[E - E(p) - EfaM) E-E(p)\ 
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Here /V 2 = + ^i),E'(p) = y/Pi + M*+y/p 2 2 + M?, E(p) = 

V~Pl + Ml + v/pf+Mf; M s = M D or M B ,M = M* D or M£, P : = 

P-— , r 2 --p-— • 

To present our results in a convenient form, which can be compared to 
experimental data, we shall use two standard variables on the Dalitz plot, 
invariant mass M nn of two pions, and cosine of the angle 9 of emitted pion 
(7r + in n + n~) q 2 = M 2 W = (u>x + uj 2 ) 2 — K 2 , and pion energies u>i,u>2 are 
expressed as 

(M + M')AM + q 2 M + M' Jq 2 - Ami J (AM) 2 - q 2 n . 

wi 2 = — =F — ttt cos 9. (89) 

4M AM q V ' 

Moreover, the term k 2 + k 2 in the first exponential in ( 1BB1) is written as 
k? + k 2 = { ((M + M')AM + g 2 ) 2 + 

ens 2 (9 1 

+ (M + M') 2 (g 2 -4m 2 )(AM) 2 -g 2 )^-- I - 2m 2 = a + 7 cos 2 9. (90) 

9 J 

Now the phase space factor looks like 



dT = d 3 k 1 d 3 k 2 5(u 1 +u 2 + y/K 2 + (MQ 2 _ _ 
(2vr) 5 2wi2cj 2 

(M + M')(M 2 + M' 2 - q 2 ) ... A 9ir „ o^i/o, , „ , , 

= 2(2M)3(2tt)3 —{[(AM) - g 2 ][g 2 - Ami}} 1 / 2 dqd cos 9. (91) 

The second exponent in (1881) can be written as 



(AM) 2 - q 2 f (M + M') 2 -q 2 \ = (AM) 2 - g 2 
4/? 2 " 4/3 2 I, 4M 2 )- 4/? 2 6) - {J2) 



Now the basic elements in (188]) (as will be shown later) which define the 
form of the tttt spectrum can be written as follows 



f d 3 p e ^I n {p)I n/ {p)p 2k = 1 {k) 

J {2-nf £! + E + AM™/ ~ (l£\ + E + AM™/ BB ' 1 J 
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' ! '' P e-^p^I n ( P )LAp)=P 2 o k Q^ 



(27T) 



(94) 



and Qnn'~ dimensionless polinomials in ratios of are given in Appendix 4. 
Substituting ([93]), flS), (l9"31) in (pi) one finally obtains. 

dw nn ,(q, cos 6) = d$\M\ 2 , (95) 

where <i$ is 

f fM br V (M 2 + M' 2 -q 2 )(M + M') r-——: r ^ ;r, 

(96) 

with 



m 2 Q + n^ 2 
2Q 2 



<EJ + ^ + AM* , 

2A/* 1 i» 



+ (f «-> 2) 



V2^ 2 y V 3 



4£» V ; 



^Ei + AM nn , 

2M nn 



+ 



+ u! + u 2 + AM nn , 



(97) 



Here we have defined ( = c n c n ,p { ®l, (tt') 3/2 , AM* n , = M B + Mb* - M(nS), 
AM nn , = 2M B - M(nS), 5 = ^-- (A ^ 2 " g2 < 1 and the parameter, which 
will be of importance, measuring the relative weight of P-wave to S"-wave 
intermediate states in the decay, 



Pr, 



rnn' 
Pnn' 



(98) 



The matrix element Ai is written for the case of the intermediate (in 
most cases virtual) channel, e.g. for T(nS) — > T(n'S)7nr this channel is BB* 
for the first square brackets, next channels are B S B* and B*B* respectively. 
In the next section we shall discuss general properties and numerical values 
of the obtained expressions. 
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7 Results 



7.1 The 7T7T spectrum 

It is convenient to rewrite dw nn > in fl95l) as the product of three characteristic 
factors: a combination of coupling constants > P nase space factor 

^ s (g + M 12 —q 2 ){M + M-) ^^^-^-^^ 

(99) 

and matrix element squared, 

dw nnl (q,cos9) = ^^ dT ph \M\ 2 . (100) 

We now evaluate all the factors in Ai, Eq. (1881) . for three typical transi- 
tions, 

A) T(2S) -> T(lS)nn 

B) T(35) -> T(lS)7T7r 

(7) T(3S) -> T(2 1 S)7T7r. 

The corresponding values of parameters (, a,j,p nn >, (p) computed with 
the SHO wave functions are given in Appendix 4. 

Before going into the details of comparison of our results and experiments, 
one can notice that the cos# dependence is relatively weak and integrating 
over cos#, and neglecting weak dependence of a, 7 on q, one can rewrite (J9j 
as follows: 

q 2 -4m 

M = a(q) - b(q) S a th - b th e~^f 
= a th - b th + b th 1 - 



^ h -b th ( q —-^)(i + q -—^ + ...) (101) 



q 2 - Am 2 \ f q 2 - Ami 

Apr) { + ~m 1 

Hence the character of the irir spectrum for small z = q is defined 

by the relation between values of A th = a th — b th , and b th where a t h = a(q = 
Im^.bth = b(q = 2m n ). 
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It is convenient also to introduce dimensionless variable x, changing in 
the interval [0,1] for all three transitions A),B),C) 



x= H - /i 2 = (AM) 2 - 4m 2 , 

fx 2 



(102) 



and n\ = °- 234 GeV2 ; /4 = °- 721 GeV2 > Mc = °- 03 Gey2 - In terms of x the 
phase space is simple 

^ = (M 2 + M- - 4^- + m i ? y ffT=T )ime ( 10 3) 



which is roughly 



jii'i 



dq 



/j, 2 y 1 'x(l — x)d cos 9. 



Finally one can rewrite the transition rate approximately and in the SHO 
basis as 



dw 



1 



dq 327TW, 2 V U J (4& 2 ) 2 



-dcos6\f x{l — x)b 



tit 



T] — X I 1 + 



8/95 



rX + ... 



(104) 

with 7] = ^ • at \™ th ■ One can see in (11041) . that the resulting itit spectrum is 
defined by the last factor, which provides three distinct types of behaviour, 
depending on the value of 77: 



1) 
2) 



dw 



dq-\J x(l — x) 



X 



dw 



2 ' 



dq-\J x(l — x) 



has a double peak and zero around x ~ 0.5. 



3) 77 < -0.5 



dw 
dq 



\J x(l — 2), no significant structures in the spectrum. 



As we shall see below, the behaviour 1), 2), 3) correspond to the transitions 
A),B),C) and this is supported by theoretical estimates of r\ in the SHO 
basis, and by experimental data. 

We turn now to the evaluation of the coefficient rj, using SHO basis, with 
a and b given in fl97|) . and averaged over cos 8, as it is done in irn spectrum. 



a = C 



m 2 Q + tf 
2Q 2 



2\ 2 



<P 2 > 



2M 



+ (1 <-> 2) 



(105) 
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Lw + AM ™' + + uj 2 + AM nn , 

(106) 

and 

(107) 

with 

C' = ( pi >, (108) 

f=<«s " >, 

exact values are defined in Appendix 4. 

Finally the form of the spectrum as given in ( 1104)1 depends on the value 
of t](ns — > n's), which can be written as follows 



„( nS _ rfS) = I 2 < e 4 " - ! (109) 

|/W(£mi' + £nn') e ^ / 

Wlth M = ^ mgiift, J ~ J ' n ) 1S g lven 111 Appendix 4. 

The values of £* n „ , £ nn / , £ nn / computed with the SHO functions are given 
in Appendix 4 and the resulting values of T)(nS — > n'S 1 ) are listed below in 
the Table 1. 
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Table 1 

The values of the parameter r], computed with SHO functions (upper line), 
with AZI restrictions (middle line) and fitted to experiment using (115) (lower 
line) . 



(nS -> n'S) 


2S — > IS 


3S ^ IS 


3S ^2S 


VSHO 


< 0.45 


< 0.27 


-3.66 


Vazi 


0.051 


0.39 


-3.2 







0.52^0.57 


-2.7 



One can see, in (105) that the magnitude of rj depends strongly on the 
value of 9, and we list in Table the maximal values of tjsho- 

Comparison with experimental values can be done fitting n to three spec- 
tra in transitions A),B),C) which yields values given in Table 1. One can 
see that the case C) agrees quite well, while in cases A) and B) the SHO 
values of rj in table 1 are in the correct ballpark, however one needs a more 
accurate calculation, since both p nn > and rj nn i depend strongly on the form of 
wave functions, indeed the overlap integrals in matrix element enter there in 
the 4th power. The obtained results and accuracy are enough only for our 
semiquantitative analysis and need to be redone with realistic wave func- 
tions. One way to obtain stable results is to impose the PCAC-Adler zero 
requirement on Ai, which is done in the next subsection. 

7.2 The PC AC improvements 

As was discussed in section 3, the cancellation between the tadpole double- 
pion amplitude (equivalent of our amplitude "b" ) and the one-pion amplitude 
(equivalent of our "a") ensures, that pion operators enter as d^ir and hence 
satisfies the Adler zero condition. In our discussion above a certain inter- 
mediate channels were chosen for a and b, different in general, and for the 
resulting approximate amplitude M, Eq. fl97|) .it is not clear whether it does 
or does not satisfy the Adler zero condition. To clarify the situation in this 
subsection we keep the form of amplitude (1971) however require the exact 
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fulfillment of the Alder zero condition. To this end we rewrite the amplitude 
(I97p in terms of momenta k 1; k 2 . One has thus the representation which will 
enable us to obtain below the "Adler-zero-improved" (AZI) form 



2 i 2 

1 



k 1+ k 2 

M = const I a e ^ I — : V 



k 1+ k 2) 2 



be \wr^ + o?> : I "I- (110) 

+ AM nn , ^L> + AM nn , + u x + u 2 1 

It is clear that when ki = u\ — 0, (M is symmetric in 1,2), indeed M 
vanishes whenever conditions are satisfied: 



1) 5(ki=u;i = 0, k 2 ) = b(ki = w x = 0, k 2 ) (111) 

2) -^U- + AM* , = ^-J- + AM nn r = r nn i. 
2M* 2M 

As a result one obtains for the tttt transition matrix element a simple repre- 
sentation (the AZI form) 

ki+k' 

M nn ' = const I e 4,3 2 I — 1 h 



(k 1+ k 2) 



4/J. 



(uJi) + LOi T nn '(uJ 2 ) + UJ 2 



+ 7— — TTTTTTT- ■ ^ 12) 



'(0) T nn '(ui + UJ 2 ) + LOi + UJ 2 



It is clear, that (11121) satisfies the Adler zero condition with any function 
T nn i(u>). We have here two model parameters: /3 2 , i.e. the radius of the 
B, B* (or D, D*) mesons, and T nn > - the average transition energy. We fix 
(3 2 by the r.m.s. radius of B meson (for T transitions) as we did before, 
(3 2 = 0.5 GeV, and r nn < calculate from the matrix element in Table 8. 
Knowing average values of given in Table 6, one immediately obtains AZI 
coefficients £ n „> = t~£ and f nn / = (r nn / + (u x + ^ 2 ) nn /) _1 . 

Finally, expressing k 2 via q 2 and COS6 1 using (I901) . (1921 one obtains the 
AZI matrix element in the form. 
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We now consider three typical cases of T(n) — > T(n')7r7r. Here the values 
of £* = t+ 1 ^^ , £ = i and £ = T+ 2^ are given in Table 8. They are computed 
as in (93), and the Adler- zero- improved t]azi can be written as 

A Ql / OC* p 2 -q — i cos 2 9 \ 

^ = ^ ( -^^e^^- - 1 ) . (114) 

The resulting 7T7t spectrum has a simple approximate form 

— — = const a/x(1 — x)|7y(nn') — x| 2 (H5) 

The spectra corresponding to (I115P with t]azi fitted to the data are 
shown in Fig. 6 together with the experimental data of CLEO Collaboration 
|42j . One can see a reasonable argument for all three types of behaviour, 
A),B), and C)with this simple parametrization. As it is seen in Table 1, the 
computed values of t/azi are very close to the fitted ones. 



7.3 Total yield of tttt 

In this subsection we derive the width of the dipion decays of heavy quarko- 
nia, T n7T (nS — > n'S) using our general equation (11001) . (11041) 



r 7r7r (nS' — > n'S) = J dw nn i (q, cos 9) 



Mbr\ 4 /I 7 2 f 1 dXy/x(l-x) |2 

1 2v b th I 7 . o \V-x\ ■ (116) 



327r 3 iv c 2 v u j m 







i 4m 2 



We shall concentrate first on the T(2S) — > T(lS)7T7r transition, since here 
the SHO wave functions might better imitate realistic ones, while for higher 
nS states the difference could be much larger. 
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Assuming first that r\ <C 1, we put it equal to zero (as it follows from the 
shape of spectrum in experiment). Then the phase space integral with x 2 
yields 

./- -y/ l-xd x =(m6; for /i 2 = 234 GeV 2 ,T(25) -> T(15). (117) 



This yields 



r-PS - is) -^7v?(x) 4 MF' ' 1166 - (U8) 



with 6 th (2S — ► IS) — 0.92C and ( is given in Table in Appendix 4, ( 2 (2S 
IS) = 0.194. As a result one obtains 



T^(2S - 15) = f^j 13 - 10 9 GeV = 0.013^^) keV. (119) 

Taking M br = 1 GeV as follows from r('0(377O) — > DD), see section 5, 
one obtains r T7T « 171 keV, while experimental value is [63] r exp (T(2S') — > 
Y(lS') 7r+7r ~) = (8.1 ± 1.2) keV. This distinction is not surprising, taking 
into account the fourth power of (j^- j and also the fourth power of overlap 

integrals between T(nS) and BB wave function. 

In particular, if one fits the value of (3\ for T(2S) to reproduce the po- 
sition of zero in the realistic x-space wave function ([181 E2], and private 
communication by the authors), then (3% will increase 1.6 times as compared 
to that fitted to the r.m.as. radius of T(2S), and the resulting value of 
( 2 and hence of T n7r (2S — > IS) decreases twenty times (!). Consequently, 
ri^ r ° f %tted ^(2S — > IS) ~ 8.5 keV, which is comparable to the experimental 
result. 

To check, whether the mass dependence of our formulas is correct, we 
now turn to the case ip{2S) — > J/tpirn. Calculating again with the SHO 
wave functions, one obtains (fitting the r.m.s. radii 0.4 fm for J/ip, 0.8 fm 
for tp(2S), and 0.58 fm for D) (3[ = 0.94 GeV, ft = Pi(ip(2S)) = 0.47 GeV, 
02 = 0.426 GeV. As a consequence (see Appendix 4 for details) one has 
pty = 2.66, and = -0.454, and b th = -0.11 GeV" 1 and the result is 

rl = 3.8 (j^-XlO- 9 GeV. 
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Fig.6 



Taking again M br = 1 GeV, as fitted to the case -?/>(3770) — > DD, see section 
5, one finally gets = 51 keV. This should be compared with PDG data, 
Tcxp = 158 keV. One can see that our theoretical value is within factor 
of 3 from the experimental result, which is reasonable, taking into account 
the crudeness of the SHO wave functions. Again, since the SHO result is 
proportional to -h ~ (t\d) 4 , where Tb is the radius of D meson, changing 

by 40% one obtains the 4 times increase of T^, making it comparable to the 
experimental value. 

We do not compare total yields for higher T(nS), with n > 3, since 
there the total SHO yield becomes too large due to specific properties of 
SHO vawefunctions, where (p n )sno grow almost factorially with n, which is 
unrealistic. 

The resulting dipionic spectra of T(nS) — > T(n'S)im with n = 2, 3 and 
n' — n — 1, n — 2, are given in Fig. 6 in comparison with the CLEO results 
from [4"2] . 

One can notice that the one-parameter fit of the form (j!12p is working 
reasonably well, yielding values of t/azi ~ ^cx P (fit) in Table 1. These values 
of r/exp are in the same ballpark as the theoretical values r] SHO given in the 
same Table 1. 

At this point it is important to stress that in the fits and estimates it was 
the value of rj averaged over angle 9 that was used, as in Eq. (I109[) . However, 
rj has the 9 dependence visualized in Eq. (11051) . and it can be written as 

,2 



-^1 = (pee ^ - 1 w (v)e + (cos 9 - cos 9)—^ 



Therefore 1) a strong 9 dependence appears when jjL is large, as in the 
3S — > 15* case, when = 0.4, and the enhancement at 9 w is explained 
seen in the CLEO data [42]. In 2S — >• IS, ~ 0.16 and enhancement is 
more shallow. 

2) The dip in the tctt spectrum in (3S — >■ IS') case is at x = r](9) and 
depends on 9. integrating the spectrum over d9 one is summing over different 
dip positions, which results in a partial filling up the dip, in agreement with 
experiment, see Fig. 6. 
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8 Discussion 



Our strategy in this paper is drastically different from the generally accepted 
"multipole expansion - Adler zero" approach. The latter is essentially ex- 
ploiting the idea of the small-size source of gluon fields. In particular the ap- 
proach of the ITEP group [TTJ EU [36] is aesthetically appealing, starting with 
the local condensate of gluonic field, turning it into the energy-momentum 
tensor of gluons and finally into that of pions, predicting in this way a peak 
at high im mass and a correct ratio of rj/im- yield in (2S — > IS). 

Our argument is that this might well be true for small size systems, such 
as toponium, however becomes inapplicable for systems of size R larger than 
the vacuum correlation length A ~ 0.2 fm [19J. Actually all decaying systems 
in question have r.m.s. radii not smaller than 0.4 fm, when nonlocality of 
gluon condensates is vitally important and brings confinement (see [13] for 
discussion). It was checked in |64j that the local gluonic condensate strongly 
overestimates binding energy of all states with R > 0.4 fm, while nonlocal 
ones (confinement) yields correct masses of T(nS) with n > 1. 

Therefore our approach is an attempt in applying nonperturbative QCD 
methods to the large size systems, using derived in this way quark-pion La- 
grangian [46]. It contains the only parameter- Mb r , which we fix by indepen- 
dent input - the width of ^(3770) -> DD, yielding M br « 2lu » 1 GeV. 

The 7T7T production in our approach is due to the presence of light qq 
pair in the body of heavy quarkonium, which should be 0(1/N C ) in am- 
plitude. Each pion is emitted with 4- amplitude, which yields nn produc- 



tion amplitude O (j^-p) = O )■ One can see in flllOl) . that indeed 



Furthermore, when the qq pair is formed inside T(nS) the latter lives part 
of time in BB or BB* (or else _B*.B*)states. The quark-pion Lagrangian (T38l) 
contains vertices with emission of any number of pions (or kaons and etc.), 
and essentially there appear two types of amplitudes: (a) with one-pion 
emission at each vertex, see Fig. 4, (b) with two-pion emission at one vertex 
and zero-pion emission at another, see Fig. 5 . 

It was shown in the paper, that due to chiral properties of quark-pion 
Lagrangian (essentially equivalent to the Adler zero requirement) these two 
amplitude enter with different signs, M = a — b. 

Finally, a and b, when computed with any reasonable (decreasing with 
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argument) functions are decreasing functions of different pionic variables, 
namely, a = a(k 2 , k 2 .), b = b((ki + k 2 ) 2 ), in case of SHO functions, a SHO ~ 

This fact implies completely different dependence in terms of q 2 = (k\ + 
/C2) 2 and cos 6, indeed, b does not depend on 9 and grows with q 2 , b SHO ~ 
exp while a SHO ~ exp ^_ cos 61 j _ The same type of behavior is ex- 

pected for realistic wave-functions, a ~ fa(oi(q), r y(q) cos 2 9),b ~ fb(q 2 ), with 
/ a decreasing and /& increasing functions of their arguments. The last feature 
is the cancellation of = a(g 2 = 4m 2 ), and 6^- the complete cancellation 
in the case A) (2S —>■ IS), a partial cancellation in case B) (35 — > IS 1 ), and 
no cancellation at all in case C) (3S — > 2S). This feature is qualitatively 
supported by the SHO calculations (see Table 1), and needs to be confirmed 
by the realistic calculation. 

As a result one obtains a simple one-parameter representation of the kit 
spectrum which works well in case of T(2S), T(3S) and if)(2S) as shown in 
this paper. Moreover, one finds explanation not only for the form of the 
spectrum, as shown in Fig.6, but also qualitatively for the cos# dependence. 

The extrapolation of the method to the cases of T(AS) decays [25], ;26j and 
to the T(5<S') transitions, found recently by Belle [65], is straightforward. One 
can easily see, that the corresponding tttt spectra are easily described by 77 
parametrization, however the SHO calculation of spectra is hardly applicable 
for such high n, n', and one needs realistic wave functions and possibly more 
channels of B, B type should be included. This work is planned for the 
future. 

9 Summary and outlook 

1. We have constructed the amplitudes, tttt spectra and the total width 

T nn for all dipionic transitions T(nS) — > T(n'S)7T7i with n = 2,3, 
and n 1 = 1, 2. The only parameter of our approach, M^. in the effective 
chiral Lagrangian, j -^(x)Mb r exp(i , y 5 (j))ifj(x)d 4: x is fixed using the decay 
^(3770) -> DD. 

2. It is shown, that dipion spectra have the form dT n7T = phase space times 

1 77 — x| 2 , where n ~ a th-hh an( j a ^ re f ers to the threshold amplitude 

®th _ 

with sequential decay of the type T(nS) — > BB*n — > T(n'S)im, and b t h 
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- the amplitude for the two-pion- vertex emission: T(nS) — > BBtttt — > 
T(n'5) plus T(n5) -> 55 -> r(riS)mr. 

3. It is shown, that 7T7t spectra appear of three kinds: A) when 1 77 j 1 the 

spectrum behaves roughly as (q 2 — 4m 2 ) 2 , as in (25 — > IS 1 ) transition in 
tp(2S) and Y(25), -B) when 77 ~ | there appears a dip in the middle of 
spectrum, as in the case Y(35 — > IS 1 ), Y(45 — > 25), C) when 77 < —0.5, 
there is a moderate enhancement at large q 2 due to amplitude b, as in 
the T(35 — > 25'). The case T(45 — > 15) seemingly belongs to case B) 
with r] cxp (AS -> 15) « 0.30 and 77^ (4S -> 25) « 0.61. 

4. We have checked the method, computing absolute normalization of 
■0(25 — > IS)itit and T(25 — > 15)7T7t using SHO functions and found 
qualitative agreement. We also computed 77(cos#) with the same SHO 
functions with parameters fitted to the known r.m.s. radii of all states. 
Remarkably, in spite of crudeness of this approximation, we have found 
a reasonable agreement with experimental Tin spectra, especially when 
the Adler zero condition is imposed (see t\azi in Table 1 compared to 

5. We have argued, that multipole expansion is not suited for all nS states 

with n > 1, since their size is larger than correlation length of gluonic 
vacuum, A = 0.2 fm. 

6. We have casted some doubt on the use of PC AC vanishing of ampli- 
tude as a sole source of damping of nn spectrum , giving example of 
amplitude with Adler zero, but actually not changing across available 
phase space. 

Moreover, the (q 2 — 4m 2 ) 2 behavior in 25 — > 15 is not due to PCAC. 
Still, relative sign of a and b is necessary for the Adler zero condition 
and is vitally important for description of irn spectra. 

7. It is argued that tttt FSI is not affecting much the form of spectrum, 
in T(nS),n < 4 since both major features: the (q 2 — 4m 2 ) 2 damping 
at threshold for (25 — > 15) and the dip in (35 — > 15) are seemingly 
not connected to FSI. However, FSI might be important for the part 
of spectrum around q ~ 1 GeV. 
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8. In developing the formalism, we have suggested the method of construct- 

ing the relativistic amplitudes of multipoint and multihadron types, the 
latter expressed via eigenf unctions of relativistic Hamiltonian. 

9. The same method is easily generalized for the case of scattering am- 
plitude of pion on heavy quarkonia, and preliminary study reveals a 
strong interaction, which may be of importance for the explanation of 
the states like Z(4430), recently discovered by Belle [66J. 

10. For n > 5 the T(5S — ► n'S) transitions proceed via the open BB 
channels, and our Eq. fl88|) predicts a strong absorptive part in the 
amplitude M. (97) which is proportional to the AtcpbMb and produces 
an enhancement factor ~ 10 3 in the total tttt width, in experiment 



r^ p (5S — > IS) ~ 1 MeV (65], whereas T%?{3S -> IS) ~ 1 keV [32]. 



This fact gives an additional support to the mechanism proposed in the 



11. As it is clear from (11191) the main dependence of total width on the 
mass of emitting Nambu-Goldstone mesons comes from the /z 7 depen- 
dence. Hence for the T(55 — > IS) transitions the ratio ^^qi^T^ ~ 



cent experimental data [65J. 
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paper. 




0.104 fa = ^(AM 2 ) -4m?), which roughly agrees with re- 
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Appendix 1 



Factorization of the Dirac bispinorial structure of the string 
breaking Green's function Gb r 

Writing in gO) J qQ = ^ Q (u)r Q $(u, v)Vq(v) and 

^ (u',v') = -^ Q (u')r Q $(u>> Q (v'), (Al.l) 

and for the light quark vertices 

Jqq(x) = $ q (x)r x ip q (x), J+(y) = $ > q (y)T 'yip q (y) (A1.2) 

one has the following structure for a given choice of two intermediate heavy- 
light bosons (?pQ(w)rQ q i{) q (w)), (ip q (w')r q Q^Q(w'), as depicted in the left 
part of Fig. 7 (the right part has a similar structure). 



w 




Fig 7 



G QQm = tr L(r Q S Q (u, w)TQ q S q (w, x)T q S q (x, w')r qQ S Q (w', v)). (A1.3) 

Here tr^ denotes the trace over Dirac bispinor indices. In the FFSR one 
writes the quark Green's function as (in Euclidean notations), k = q,Q 

S k {x,y) = {mk-b^ml-D 2 )- 1 = 

poo ps 

= (m k -D) ds(Dz) xy e~ K $(x,y)exp(g dra^F^) (Al.4) 
Jo Jo 
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The last factor in (IA1.4I) takes into account spin-dependent interaction 
and can be disregarded in the first approximation. Then all Dirac matrix 
structure of (IA1.4I) is furnished by the prefactor (m k — D), which can be 
rewritten in terms of momentum as shown in [57] , namely for a simple quark- 
antiquark loop one has 

Y T = ^tr L (m 1 -JDi)r(m 2 - A2); r = ^tr L ((m 1 -ipi)T(m 2 + ^ 2 )T. (A1.5) 

Here p t = p^7^, i = 1,2, and in the cm. system p^ = —pf^ = Pi, = 
iu>i, = iu>2, so that 

Yr = ^tr L (T(m 1 + ^174 - ip k ^ k )T{m 2 - ^74 - ipiji))- (Al.6) 

Note, that the index 1 is referred to a quark, while index 2 to an anti- 
quark. In (IA1.3I) one has heavy quark Q with momentum p^ = (iujQ,p^ ), 

heavy antiquark Q with momentum p®' = (iu}Q,p^), light quark q with 
momentum = (io> g , and light quark q with momentum q^ = {iuj q ,qi). 
Writing S k (x,y) = (m k — D)G k (x,y) and taking all G k out of the sign of 
trL, since they are proportional to the unit matrix in Dirac indices, one has 
G QQqq = ZIl kG k , where 

Z = tr L [(T Q (m Q + ^q7 4 - ipf' Qq {m q - ^74 + iq k ^ k ) 

T q (m q + uj q ^ 4 - iqai)T Qq (m Q - ujq^ a + ip^ln))]- ( A1 -7) 

Matrices Tj, i — Q, Qq, q, qQ, define the quantum numbers of participat- 
ing particles, while T q can be only of two types: T q = 1 for a qq vertex without 
NG emission or emission of even numbers of NG meson and r„ = ^75%^- for 

one NG meson emission, and — ^ aX ^ bXb for two NG meson emission. 

To normalize the factor Z properly, one should take into account, that the 
path integral for the scalar part of the quark Green's function G k (x, y), k = 
q, Q can be rewritten as follows [57] 

G k (x,y)= / ds{D A z) xy e- K §{x,y)= / V JXy Due- K <f>(x, y), (A1.8) 

JO J ZuJ k 

where u k is an averaged energy of the quark, uj k = (yrn|~+p|),and for the 
white quark-antiquark Green's function one has 

■^(tr c G k {y,x)<5>{x,x)G k {x,y)${y,y)) = — — (x,x\e~ HT \y, y) = 
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J2^n(x,x)^t(y,y)e- EnT . (A1.9) 



Therefore it is convenient to introduce the projection factors where 
subscripts (+) and (— ) stand for quarks and antiquarks respectively 

(k) 

A* = 7L , k = g ,Q, (ALIO) 

zu>k 

and the normalized factor Z looks like 

Z = tr L (T Q A+T Q ,A-T q A^T qQ AQ) (Al.ll) 

Using (IM.llj) and (IA2.23P arrives at Eq. (ggj 

As a next step one separates c.m.momenta and introduces in the QQ 
system Green's function with zero cm. momentum and in the energy rep- 
resentation G^q°(t, r', E) as in (T4"6l) and the same in the QQ system. Note 

that dimension of Gqq is [Gqq°(t, r', E)] = m 2 . 

Separation of cm. and relative coordinates goes as follows. Consider 
QQqq system represented as a system of two bound states Qq and Qq with 
cm. coordinate Ri and R2 and cm. momenta Pi and P2. One can find total 
cm. and relative momentum of the QQqq system as P = Pi + P2, p = 
I (Pi — P 2 ), and cm. and relative coordinates of equal mass Qq and Qq 

T3 I > 

mesons as R = ^ 2 , r = R x — R 2 , so that one has P1R1 + P2R2 = 
PR+pr. 

On the other hand, cm. coordinates are found via the average energies 
of quarks, so that in the initial and final states of the amplitude depicted in 
Fig. 7 one has 

UJl+Ul UJ 2 + "2 

r; = ^±§^, r' 2 = ^# (A1.13) 

UJi + Hi UJ2 + H2 

For the same masses in Qq and Qq, u% — u 2 , fli = ^2 one finds 

„ u;x fiu + v ft , . , , , K , . 

R =— 75 + — 75" 5-. r=— ^(u-v =c(u-v A1.14 
to + il 00 + il 2 00 + il 

As a result the total integration element is 

d 3 xrf 3 urf 3 v = (pp) d 3 Rd 3 rd{x - u) = 
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d 3 Rd 3 (v-u)d 3 (x-u). (A1.15) 

The factor Z and 2/123 in Eq. (j4"5|) contain effective energies u and Q. We 
now explain, how these can be computed, solving spin-independent Hamil- 
tonian. 

We start with the Qq system in its cm. system. One has (neglecting 
spin-dependent forces) 

Hxp = (y/q 2 + m\ + ^/q 2 + m% + A + V(r qQ ))4> = M Q $. (A1.16) 

Here A = A SE + A string , V(r q Q) = ar q Q + V cou{ and u = (^q 2 + m 2 Q ), 
Q ee ( ^q 2 + m 2 ) . Another way is to use the einbein representation [661167] 

/ q 2 + m 2 u q 2 + m% Q \ 

(^7^ + 2 + + 2 + A + F(r ^ } J ^ = MQg " (fi ' (AL17) 

with the subsequent optimization of MQq{Vt, uo) using equation ^ = ^ = 0, 
which yields Q = Qq, uj = ujq and it was shown that ujq, Qq are close to the 
previous definition in terms of averages [67], see [68] for a review of string 
Hamiltonian technic. 

The spin-dependent forces are treated also in the FCM, as shown in [68J. 

The resulting spectrum for QQ systems was obtained in [62J, while that 
for Qq systems in [57J. It is important, that the method contains minimal 
number of parameters: current (pole) quark masses, (m c = 1.4 GeV, mb = 
4.8 GeV, m u = m d m 0, m s = 0.2 GeV) string tension a = 0.18 GeV 2 , and 
a s (q) = with M B = 1 GeV and A§££ = 0.22 GeV In some recent 

K QCD 

applications also the flattening of the confining potential due to light pair 
creation was taken into account [69], but we shall neglect this effect here. 

With the given input parameters the resulting spectrum and Q, uj are 
given in Table 1 for charmonium, Table 2 for bottomonium and Table 3 for 
some heavy light mesons. Please note, that values of Q,u depend on the 
system, where a given quark (antiquark) enter. 
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Table 2 



Values of f2 and (p 2 ) for 1 states*) of charmonium, computed in FCM in 
|18[ EH] together with masses vs experimental data 



State 




1>(2S) 


^(3770) 


V>(3S) 




Q R , GeV 


1.58 


1.647 


1640 


1.711 


1.770 


(p 2 ) R GeV 2 


0.569 


0.820 


0.746 


1.064 


1.146 
1.300 


mass, MeV 


3.090 


3.675 


3.800 


4.094 


4.442 


theory M cog 












mass, MeV 


3067.8±0.3 


3674.1±1.0 


3771.1 ±2.4 


4039 ± 1 


4421 ± 4 PDG 


exper M cog 






6 = 105° 




4411 (Belle) 
4361?( (Belle) 


( P 2 )ea GeV 2 


0.522 


0730 


0.731 


0.940 


1.141 


n EA , GeV 


1.60 


1.65 


1.65 


1.70 


1.75 



*>Qr, (p 2 )r are computed via Salpeter equation, Eq. (1A1.16I) while £Iea, (p 2 ) ea 
are computed in the einbein approximation, Eq. (lA1.17j) . 

Table 3 



The same as in Table 1, but for 1 states of bottomonium, as calculated 
in [TBI ESI- rrib = 4.823, M B = 0.95 A 5 (MS) = 0.220; (all in GeV) 

a = 0.178 GeV 2 



State 


T(1S) 


T(25) 


T{3S) 


T(45) 


T(5S) 


Q, GeV 


5.021 


5.026 


5.056 


5.088 


5.120 


(p 2 ) GeV 2 EA 


1.954 


2.026 


2.334 


2.674 


3.012 


mass, GeV 


9.460 


10.010 


10.356 


10.633 


10.873 


theory c.o.g. 












n 3 Si mass, GeV 


9.460 


10.023 


10.355(1) 


10.579(1) 


10.865(8) 
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Table 4 



Values of u>, Q and (p 2 ) for heavy-light mesons (masses computed in [57] 
also given vs experimental data from PDG) 



Meson 


D 


D s 


B 


B s 


u, MeV 


507 


559 


587 


639 


tt, MeV 


1509 


1515 


4827 


4830 


(p 2 ), GeV 2 


0.273 


0.290 


0.359 


0.383 


Mass, MeV 


1869 


1967 


5279 


5362 


(theory) 










exper., MeV 


1869.3 


1968.2 


5279.0 


5367.7 




±0.4 


±0.4 


±).5 


±1.8 



Appendix 2 

Normalization of decay amplitudes in the relativistic represen- 
tation 

We start with the amplitude of two-body decay where initial and final 
states are created by the current operators: initial state by j\ = ipTiip and 
final states ji = ipTiip, i = 2, 3, while the pair creating Lagrangian L = 
M br (ipT x ip) acts at the point x. 

The Green's function for this amplitude in the space-time can be written 

as 

3 

Gi 23 , = (0| \{m = tr(r 1 S(l,2)r 2 S(2,x)T x S(x,3)r 3 S(3,l)). (A2.1) 
i=i 

Here trace is over color, flavor and Dirac indices, and S(i, k) is the quark 
propagator from the point i to the point k, see Fig. 3. It is easy to see that 
dimension of J d^xGyizx is m 9 . To write propagators in convenient relativistic 
form, we shall use the fact, that between all quark lines in the diagram in 
Fig. 3, corresponding to the amplitude flA2.ip . acts confinement, and therefore 
all propagators are at the boundary of the film with string tension a, which 
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creates effective mass u q for all quarks, including massless ones. This fact 
was established in [7U], where this mass was shown to be what is called the 
constituent mass, and uj q was computed repeatedly and accuracy was checked 
in [SZ]. 

The quark Green's function in the Fock-Feynman-Schwinger representa- 
tion has the form (in the Euclidean space-time) 

S(x,y) = (m-D) x G(x,y), (A2.2) 



poo 

G(x,y) = / ds(D 4 z) xy e- K $ a (x,y), (A2.3) 
Jo 

as it was shown in [57] , one can rewrite G(x,y) identically as follows 



G{x,y) 



Dlv{D 3 



2uj 



e K $ a (x,y) 



(A2.4) 



and the functional integral (Du) has the meaning of the averaging procedure 
over all possible values of u = l^p, where z^(t) is the (Euclidean) time on 
the quark path at the proper time r. Note that kinetic energy is 



K 



— ) dr + m 2 s 
dr 



oo 
2 



dz, \ \ in- 
dz 4 J I 2uj 



dt 



The averaging Du with the weight exp(— K)<f> a finally yields uj - constituent 
mass. 

Therefore S(x, y) can be written as 



S(x,y) = {Zg(x,y)) U} , 



where Z 



(rn-D) 
2ui ' 



and 



g(x,y) = / (D 3 z) xy e K ^ CT (x,y). 



(A2.5) 



(A2.6) 



Consider now the qq white system, e.g. the current-current correlator 

(j(x 1 )j 2 (x 2 )) = (trT 1 S(x 1 ,x 2 )T 2 S(x 2 ,x 1 )). (A2.7) 



47 



One can write, (neglecting spin splittings and hence second exponent in 
(jih) = (Y J (D 3 z) XlX2 (D 3 z) XUX2 e- K - R (W(x 1 ,x 2 )) A ) UJt . (A2.8) 

Here 

- = trTijjn t - D 1 )Y 2 {m 2 - D 2 ) _ 



2ui2uj 2 

1 



tr(T 1 (m 1 - ip\)T 2 {m 2 + ip 2 ))- (A2.9) 
/ ko\U 2 

and (W(xi, x 2 )) is the Wilson loop average, with the closed loop along tra- 
jectories of quark and antiquark. 

In (IA2.8[) one can introduce relative and cm. coordinates r), p, and finally 
write 



(D d z) XlX2 (D- i z) XlX2 e- K - K (W) A = J (D 6 r ] ) 0fi (D i p) XuX2 e- K - K (W) A = 

= J i!L e *P(xi-x 2)(0 | e -H(x 24 -x 14) | 0) = 

= / ^r 3 ^ iM2) £ l^(o)|V^ AT (A2.10) 

Our important task is to go over from the point-point correlator like 
O1J2), to the hadron-hadron amplitude, which can be done introducing the 
so-called decay constant /p. For the standard normalization of the hadron 
state (one hadron in the volume V3 = 1) one has^l 



p - 0> - w- (A2 - n) 

with By = 1 for scalars and e-p = £^ for vectors. 

Hence one can rewrite (IA2.8D using (IA2.10D . (IA2.11D as 

U1J2) = YJ2J |^e iP(Xl - X2) - E " AT £ is 2 ^^. (A2.12) 



2 Note, that normalization of state \n) is (|A2. 1 1|) corresponds to (|A2.10j) and differs 
from b7\. 
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Comparing (IA2. lOj) and (IA2. 12j) one obtains /f - , (see [FT] for details and 
numerical estimates) 

(/ „ )2 = 2JupW (A , 13) 



n 



In a similar way Eq.( 1A2.lTl) will help us now to define hadron-hadron 



amplitudes by amputating the "current-at-a-point" matrix elements. 

To this end we represent Gi23 X as follows (we omit integral signs for 
brevity) 

GW = Z 123x (D 3 z) 12 (D 3 z) 2x (D 3 z) x3 (D 3 z) 31 \\Duj lk e" K ^{W) (A2.14) 

Here the subscript {ik) runs over (12), (2x), (x3) and (31), and we shall 
omit the path integrals over Du>i k and replace it by the sign of averaging 
over uj of the whole expression, G, e.g. (G) u , which leads to the appearing 
of average values u in Z, and the (multiparticle in general) Hamiltonian H. 

To introduce this Hamiltonian and corresponding eigenfunctions, one 
should define the hyperplane, and we choose it as a hyperplane at the time 
point x 4 , which crosses path (12) at space point u and path (13) at the point 
v. One can divide the paths at these points: 

{D A z) 12 = (D 3 z) lu d 3 u(D 3 z) u2 , (D 3 z) 13 = (D 3 z) lv d 3 v(D 3 z) 03 (A2.15) 

and introduce cm. and relative coordinates in three regions (on three pieces 
of the Wilson plane), denoted in Fig.3 by letters a,b, and c. For the cm. inte- 
gral (-D 3 p)y,z one can use as in (1A2.10P the representation J exp(iP(y — 
*))• 

Denoting the cm. and relative coordinates in regions a,b,c as X±, X 2 , X 3 
and r]i,r] 2 ,r] 3 respectively, one can rewrite Gi 23x as 

G l23x = {Z U3 ^D 3 Vlk ^d 3 ud 3 ve^^-^x 

x e -iP 2 (X 2 -X ra )-*P 2 (X 3 -X^ e -^^}^ (A2.16) 

At this point one can use as in (IA2.10I) the connection of the einbein 
Hamiltonian with the path integral (see [13] and paS] for details) 

J (D% k )r u r 2 e- K *(W ik ) = ( ri | e -^ T |r 2 ). (A2.17) 
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Here and defined for a pair of quark paths i and k, are 



F ife = + ^ + + ^ + V%*, L 2 ). (A2.19) 

2 2^ 2u; fe 2^^ V ' ' ; V ; 



Here un- = ,and the string potential V was derived in [70] and given 
in Appendix D of [57] . 

The r.h.s. of ( 1A2.17I) can be written as a spectral sum 

( ri |e-^ T |r 2 ) = J2Mn)^n(r2)e- EnT . (A2.20) 

n 

Using (1A2.16I) and keeping only fixed states n!,n 2 ,n 3 in regions a,b,c 
respectively, one gets for the Fourier transform of Gi 2 3 X 



Gi 23 ^ 4 a;)p ii p 2i p 3 = Z 123x J c/ 4 xc/ 3 Mrf 3 t;^ ni (0)^ i (u-v)^ n2 (u-x)^; 2 (0) 

x^ n3 (x - v )^ 3 (o) e - iPlX - +iP2X - + ' jP3X -x 
exp{-E ni (t! - Q - E n2 (t x - t 2 ) - E n3 (t x - t 3 )}. (A2.21) 



x 



With the definition of cm. coordinates R, X w etc. as in (1A1.14[) . one can 
write 

d 4 xd 3 ud 3 v = d 3 Rd 3 (v - u)d 3 (x - u)dt x (A2.22) 

integrating over d 3 v one gets (27r) 3 5^ 3 )(Pi — P 2 — P3), and integrating over 
dt x yields the factor 2tt5(Ei — E 2 — E3). 

Now we are in position to amputate the "current-at-a-point" pieces ip ni (0) 
and go over from point-to-point correlators (Green's functions) to the hadron- 
hadron amplitudes. To this end one must replace every current vertex 
(0|j;A...|0) by (0|j i |n i Pj)(n i P i |A...|0) and delete the factor on the left. More- 
over, one deletes the energy factors exp(— £\ E n .ti). As a result one obtains 
hadron-hadron amplitude 

(n 1 P 1 |G|n 2 P2,ri 3 P3) = 
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(/G , i 23a; ^)p li p 2i p 



3 



(2tt) 4 



where we have defined 



5 (4) (Pl ~V % - Vz)J ni , n2> nMn 3 (A2.23) 



1/123 = — r =^=, (A2.24) 



2,3 ~ 4 



Jnm 2 n 3 (P) = 1/123 ^ ^ (v - u)rf 3 (x - U^P'^m (u - v)^n 2 (u - x)^„ 3 (x - v) . 

(A2.25) 

Here p = |(Pi — P2), r = c(u — v), and c is defined in (1A1. 14[) . The 

values of Zs for different Tj = 5.V, A, P can be obtained as Z = Yri , and 
were calculated in [57] . 

Finally, the decay probability can be written as 

iw J^ {Vl . V2 . nVm ^^. (A2.26) 

One can check, that dw and r = J dw have correct dimension, if taking 
into account that dim (1/123) = dimT x = [m], and dim(J nin2n3 ) = [m~ 3 ' 2 ]. 

We note in conclusion, that one could define another normalization of the 
bound states \n, P >, such that 

d 3 P 

AB = > Aln. P >< n, PIS = > A\n >< n\B. 



In this case the new amplitude 

(TiiAlGl^naA) = x/2^Z^2^(n 1 P 1 |G|n 2 P 2) n 3 P3) 

and in the definition of probability one should divide (1A2.26|) by the factor 
f| 3 =1 (2_E n J, which amounts to same dw and V. 

It is easy to see, that the emission of any number of additional pions leads 
to an additional dimensionless factor IT. cx p^' x = j n the amplitude J n , non , 

and appearance of additional dimensionless factor Yli ^fp' i n dw, therefore 
all normalization stays intact. 
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Finally we can compare dw in (1A2.26I) with our expression (47) for T n . 
Identifying J nin2m with ( v /7«/nn 2 ,n 3 (p)) in (!4"81) . one can see that normaliza- 
tion of both expressions coincides, one can also check, that in the nonrel- 
ativistic limit 3/123 = ^fyZnk, however for decay of heavy quarkonia into 2 

heavy-light mesons, 3/123 is twice as big in the static limit due to a/ Z 2 Z 3 in 
( IA2~24D . 



Appendix 3 

Kinematics of the dipion decays 

We start with the two particle intermediate state, eg. BB or BB*, and 
define momenta of each hadron 1 or 2 as Pi, P2 and the total cm. momentum 
is P = Pi + P 2 . Then one can write Pi = Pq + p q , P 2 = Pq + p„, and 
introduce relative momenta in hadrons 1 and 2: q 1 = |(Pq — Pg), q 2 = 

!(Pq-p 9 )- 

In the total cm. system one has Pi = — P2 = p, and hence p q = |p — q\, 
p q = — |p — q 2 . Therefore the combination entering Z in (1551) . p q — p q = 
— p+q 1 — q 2 . This is used in f[5T?j) and one also obtains in (155]) that q x — q 2 = 0. 

Now consider the hypersurface at time t x , when the qq pair is created. 
One has quarks Q at coordinate u, qq at x and Q at v. The tctc system is 
also created at x with total momentum K = ki + k2, and the total cm. 
momentum P = P2 + K and we keep notation p for the relative momentum, 
in the BB system, p = |(Pi — P2)- We define the cm. coordinates of 
hadrons Ri,R,2 and the total cm. coordinate R, and relative coordinates 
r = Ri — R 2 , p = x — R, where 

_ ljix + ^iu _ u 2 x + f^v R _ giRi + E 2 R 2 + (1^(1) + ^(2))x 

l " wi+fii ' 2 ~ ~ u 2 + VL 2 ' ~ E 1 + E 2 +u*(l) + co w (2) 

(A3.1) 

Here U\ = ujg,u 2 = uj q and similarly for Qi,Q 2 . Writing the exponent of 
plane wave free motion of hadrons 1,2 and two pions in terms of cm. and 
relative momenta p, and R, 

P1R1 + P 2 R 2 + Kx = PR + pr + kp (A3.2) 



52 



one arrives at the expressions 

a 2 Pi-aiP 2 1 Ui + VLi . 
P = ; , k=- K, oa= — ,2 = 1,2 (A3.3) 

and a 3 = 1— «i— a 2 , where is the total energy, = E 1 +E 2 +uj 7T (l)+u) n (2). 
In the equal mass (and energy) case one arrives at the same equations as 
before, a x = a 2 , p = |(Pi - P 2 ) and p g - p 5 = q : - q 2 - p. 

However now integrating as in (jHJ) over the coordinate x where the qq and 
mi pairs are emitted one has instead of previous case without pion emission, 
Qi — = K(we neglect small corrections of the order of 

Hence for the diagram with the 2n emission one has Z ~ (p q — p q )i = 
Ki — pi as used in (1771) . 

In the second part of the diagram (the r.h.s. part of Fig.) pions are not 
emitted at the point y, and q' x — q' 2 = 0, hence there Z is the same as in the 
pionless decay case. 

We turn now to the 3 body phase space of X{n')7i7i and useful coordinates 
on the Dalitz plot. We choose as such the standard quantities: invariant 
dipion mass = q 2 = (k\ + k 2 ) 2 = (loi + 0J2) 2 — K 2 , and cos#, where 9 is 
the angle of tt + in the cm. of 7i + 7t~ with respect to the direction of incident 
quarkonium. One has 

n __Q M' 2 — M 2 — g 2 + AMu 2 

cosO = — = — ; A3. 4 

^ q 2 _ 4m 2 ^ M 2 Z. ( g + M') 2 ][M 2 - (M' - q) 2 } 

The pion energies uji and u 2 can be written as 

~ „ (M + M')AM + q 2 (M + M') Jq 2 - Ami J (AM) 2 - q 2 

vi 2 = cTd cos = — =p — cos 9. 

AM AM q 



Therefore the sum toi + to 2 and K 2 are 



(A3.5) 



q 2 + M 2 — M' 2 2 ((AM) 2 - q 2 )((M + M') 2 - q 2 ) 
Ul+U2 = 2M ' = AM 2 (A3 - 6) 



and the combination appearing in 

k 2 + k 2 = 2(£ 2 + d 2 cos 9) - 2m 2 = a + 7 cos 2 9 (A3. 7) 

where we have defined for the reaction X[n) — > X(n')7T7i; M = mass of X(n), 
M' = mass of X(n'), AM = M — M' . 
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Appendix 4 

(k) (k) 

The overlap integrals J^ n , and p^, 

One starts with the SHO wave functions, which can be written as ^f n {q, 13) = 
V n (q)e- q2 / 2/3 \n = 1, 2, 3, .. for (nS) states 

^■(^(T-»(fi ,+ S)V-^- 

The overlap integrals of the QQnS state and the n — 1 state of heavy- 
light mesons (BB or BB* for T(nS')), are 

c 2p2 /■ ^3„ 

e" — / n (p) = J -~^^ n (q + cp,^l(q,/3 2 ) = 

{}. = i;{}. = »-|^; « = ii^»l (A«) 



ft / V ft 

According to (|93l) the overlap integrals J 7 ^ with n, nf = 1,2,3,4 for 
1S,2S,3S,4S states respectively can be written as 



= / ^~^»(p)^(p) 



= c n cAtt') 3/2 P 2 o k p { nl» = c 2 = ^, c 3 = ^ (AAA) 
Using I n (p) defined in (!85|) . one obtains 
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P ( S = -|(-32/V + 6tV - 3y't* + 3yH 2 - lOytH' 2 + 7t'H% (A4.5) 

pS = -\{-v + In p5? = f (7t 2 - 3,)(t 2 - y) 

p£> = -^(-32/V + lOyy'f + 7t 4 + 5yY 2 + ^tY 2 + 21fV) 

Here we have defined in p$: p x = P(3S);p[ = /3(2S), f3 2 = /3(B) « 
P(B*) or £>,£>*. 



2/? 2 - /5 2 , 2# 2 - /3 2 



2/ = A > 1/ 



/? " 2 = A- 1 + A'- 1 



A A' 

For pij Pi = Pi@S), P[ = Pi(lS) etc. 

We are defining all Px(nS) through m.s.r. of the corresponding states, 
using SHO wave functions, fti(nS) = y ^ip ' anc ^ taking ( r2 )nS f° r T from 
HHIEB] one has (r 2 ) n5 = (0.2 fm) 2 , (0.5 fm) 2 , (0.7 fm) 2 , (0.9 fm) 2 , (1.1 fm) 2 
for n = 1,2,3,4,5 respectively, which gives Pi(nS) = 1.22 GeV, 0.75 GeV, 
0.67 GeV,0.61 GeV, 0.56 GeV. 

For B, B*(n = 1) one has (r 2 )i S = (0.5 fm) 2 , [57] and therefore (3 2 = 0.49 
GeV; for D,D*((r 2 ) D = (0.58 fm) 2 and p 2 = 0.42 GeV. 

Po is defined by both nS and n'S, therefore Po(3S, IS) = 0.92 GeV, 
Po(2S, IS) = 0.96 GeV, p (3S,2S) = 0.79 GeV. 

p m 

Now one can define the parameter p nn i = -Sf for n,n' = 2,1; 3,1; 3, 2. 

p i 

nn 

Using (1A4.5I) and values of $ quoted above, one has 



_3 (-y + ft 2 ) _ _(M±Zf)-533 
921 ~ 2 (-y + t 2 ) " 3 ' 74 ' P31 -2(- 2/ + ^)- 5 - 33 

p 32 = 6.79. (A4.6) 

To estimate fl88|) one extracts AM* n , and AM nn / given in the Table 1 for 
the T(nS) — > T(n'S)inr transitions with 55* and BB intermediate sates 
respectively. 
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Table 5 

The values of AM* n , and AM nn /(in GeV) for lowest nS states. 



State 


2S 


3S 


4S 




0.58 


0.25 


0.025 


AM nn i(E ns ) 


0.54 


0.205 


-0.02 



Note, that another nearest threshold for AM nn / is B*B*, which increases 
AM nn i by 0.09 GeV. For a rough estimate of (coi), (u 1 + u 2 ) we can use the 
values of anailable phase space AM = M n — M n > for 2S — > IS, 3S — > IS and 
3S — > 2S* cases respectively, and the values of (u^) are given in the Table 6. 

Table 6 



decay 


2S -> 15 


35^ 15 


3S ^2S 


AM, GeV 


0.56 


0.895 


0.332 


(wi), GeV 


0.28 


0.44 


0.166 




0.96 


0.84 


0.626 


a 


0.12 


0.36 


0.015 


7 


0.16 


0.4 


0.055 


{c n C n 'Pnn) 2 


0.39 


0.257 


0.17 


Pnn' 


3.74 


5.33 


6.79 



Table 7 

Parameters of the T(nS) SHO eingefunctions and overlap integrals fitted to 
the known r.m.s. radii of nS states. 



nS 


is 


2S 


3S 


A 


1.22 


0.75 


0.67 


A 


3.22 


1.375 


1.15 


y 


0.384 


0.636 


0.56 



Table 8 

The same as in Table 7 for intermediate quantities and the final one 77 in two 
lines; t/sho f° r direct calculation, Eq. (108) and r/g^Q the AZI -improved 
values, as discussed in the text. 
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nS — > n'S 


21 


31 


32 


t 


1.07 


1.08 


0.92 


Pa 


0.96 


0.84 


0.626 


(tt>y 


0.5 


0.416 


0.5 


c 


0.2 


0.107 


0.085 


r 


0.638 


0.657 


0.813 




0.824 


0.969 


1.033 




0.56 


0.5 


0.761 


VSHO 


< 0.45 


< 0.27 


-3.66 


'ISHO 


0.051 


0.39 


-3.2 







0.52 


-2.7 



Figure captions 

Fig.l. Light quark pair creation inside heavy quarkonia. 

Fig.2. Emission of two pions from the light quark loop: a) two-pion emission 
b) successive one-pion emission. 

Fig.3. The connected 4-point Green's function for the decay of heavy quarko- 
nium into two heavy-light mesons. 

Fig.4. Pictoral image of the successive pion emission with heavy-light mesons 
in the intermediate state. 

Fig.5. The same as in Fig. 4, but with two-pion emission at one point. 

Fig.6. Experimental data of the Cleo Collaboration from [42] with the the- 
oretical parametrization as in Eq. (115), ^ = const (phase space) 
\i] — x\ 2 , with n(fit) given in Table 1 in comparison to theoretical pre- 
dictions. 

Fig. 7. A diagram of two heavy-light meson contribution to the heavy quarko- 
nium Green's function. 
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